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Wetting and contact-line effects for spherical and cylindrical droplets on graphene layers: A 

comparative molecular-dynamics investigation 

Giulio ScocchiQ Danilo Sergi, Claudio D'Angelo, and Alberto Ortona 

University of Applied Sciences (SUPSI), The iCIMSI Research Institute, Galleria 2, CH-6928 Manno, Switzerland 

(Dated: November 22, 2011) 

In Molecular Dynamics (MD) simulations, interactions between water molecules and graphitic surfaces are 
often modeled as a simple Lennard- Jones potential between oxygen and carbon atoms. A possible method 
for tuning this parameter consists of simulating a water nanodroplet on a flat graphitic surface, measuring the 
equilibrium contact angle, extrapolating it to the limit of a macroscopic droplet and finally matching this quantity 
to experimental results. Considering recent evidence demonstrating that the contact angle of water on a graphitic 
plane is much higher than what was previously reported, we estimate the oxygen-carbon interaction for the 
recent SPC/Fw water model. Results indicate a value of about 0.2 kj/mol, much lower than previous estimations. 
We then perform simulations of cylindrical water filaments on graphitic surfaces, in order to compare and 
correlate contact angles resulting from these two different systems. Results suggest that modified Young's 
equation does not describe the relation between contact angle and drop size in the case of extremely small 
systems and that contributions different from the one deriving from contact line tension should be taken into 
account. 
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I. INTRODUCTION 

Since the discovery of buckminsterfullerenes yl (awarded 
with the 1996 Nobel Prize in Chemistry) and carbon nan- 
otubes [2], carbon nanoparticles have been the subject of in- 
creasing scientific interest. In 2010, another Nobel Prize (this 
time in Physics) was assigned to Geim and Novoselov JH] 
for their research on graphene, the basic building block for 
graphitic materials |@], where graphitic materials are defined 
as any allotropic form of carbon characterized by a six atom 
ring structure. Attention to this kind of materials is justified by 
their peculiar properties as well as their possible applications, 
which range from electronics [5] to polymer-based compos- 
ite materials [6]. Amongst the different challenges presented 
by carbon nanoparticles processing, one of the most demand- 
ing consists of achieving fine dispersions of these materials; 
this subject has indeed been thoroughly investigated in recent 
yea rs, in experimental 171- 11 111 as well as in modeling studies 
I12 - [lql . In simulation studies involving aqueous dispersions 
of graphitic materials, modeling the interactions between wa- 
ter molecules and carbon atoms obviously represents one of 
the most important tasks. These interactions, regardless of 
the water model, are usually represented as a simple Lennard- 
Jones pote ntial between water oxygen and graphitic carbon 
atoms lIUll: 



V(r)=4eco 



\ to ' 



v to ' 



(1) 



where Geo determines the equilibrium distance, £co is the 
depth of the potential well and rco is the distance between 
a pair of oxygen and carbon atoms. Non-bond interaction 
between hydrogen and carbon are usually not accounted for, 



even if some authors have investigated this option 112411 . One 
way to proceed to estimate suitable values for these parame- 
ters consists of simulating a water droplet on a graphitic sur- 
face and comparing the resulting equilibrium contact angle 
to experimental measurements, it being understood that the 
modified version of Young law for wetting has to be 

considered in order not to neglect the influence of contact line 
tension at the nanoscale 1331: 



cos 9 = cos doc 



(2) 
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In the above formula, 9 and are the average actual and 
macroscopic contact angle respectively, K is the contact line 
tension, y is the superficial tension of water and r ca is the con- 
tact area radius. In any case, it has to be noted that wetting 
at the nanoscale still remains a largely unanswered question 
Il32ll . to the point that it is not trivial to even define the con- 
cept of contact angle in this size range [34]. A similar un- 
certainty affects the definition and estimation of contact line 
tension, whose influence on contact angle is considered sig- 
nificant only for droplets with a diameter ranging from one 
to some hundred nanometers, according to different authors 
lHjHHI]. Even the sign of line tension is not unambiguous and 
some authors even hypothesized it to be a mere artifact due 
to poor experimental measurements, which are indeed known 
to be quite complex and sensitive [32fl. For a more detailed 
description of these issues, we refer to the excellent reviews 
by Mendez- Vilas et al. H32I1 and Amirfazli and Neumann 13511 . 

In molecular dynamics studies, estimates of k were ob- 
tained by simulating nanodroplets of different sizes, plotting 
contact angles as a function of droplet base radii (for a fixed 
set of Lennard-Jones parameters) and extending the relative 
linear fits for l/r ca — > 0. This approach was followed by 
Koumoutsakos and co-workers in an accurate series of papers 
lfl8[ I24I l25ll . which constitutes a reference for MD investi- 
gations of the same nature. In their first study, along with a 
comprehensive review of water-carbon interaction potentials 
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for different water models, the authors provided a detailed de- 
scription of the procedure (originally developed by Blake et 
al. IDal and de Ruijter et al. IB7I0 used for retrieving the av- 
erage contact angle from MD simulation of droplets on sur- 
faces. Moreover, they found a direct linear relation between 
contact angle and water monomer binding energy on graphite 
(actually modeled as two graphene sheets) and determined 
the values for £co and afco which matched the experimen- 
tal contact angle of water on graphite 01811 . In the following 
work, they explored the influence of Lennard- Jones cutoffs 
and hydrogen-carbon interactions on contact angle and com- 
pared a discrete graphitic plane model with a continuous one 
Il24ll . In the last paper of the series, they investigated the influ- 
ence of fluid and surface impurities on the wetting properties 
of water on graphite [25]. 

The values for Eqo and afco proposed in these papers 
Ell HU, l25ll have been used for similar MD studies by many 
researchers OEMIlIrlS- These two parameters have been 
adjusted in order to reproduce contact angles of water on 
graphite ranging from 42° to 86°, according to what was re- 
ported in the literature [|38Tl4lll and generally accepted by the 
scientific community [42[| . However, these experimental mea- 
surements are at least thirty years dated. Importantly, they 
were performed on graphite samples whose surface character- 
istics and purity, even in the case of pyrolytic graphite ll39tl 
and highly-oriented pyrolytic graphite 11401 14111 . can hardly be 
compared to the ideal atomic smooth graphene sheets used 
in atomistic simulations to model graphitic surfaces j43ll . In 
2009, an experimental measurement closer to the ideal con- 
ditions of the simulations reported above was performed by 
Wang et al. l44ll . who produced samples consisting of few 
graphene layers and characterized them with WXRD, AFM 
and FTIR methods. In this study, the contact angle of water 
on graphene layers turned out to be 127°. 

The measurements performed by Wang et al. are the most 
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Figure 1: Schematic drawing of spherical droplets and cylindrical 
filaments in the initial configuration, (a) and (b), and after equilib- 
rium is reached, (c) and (d). 



recent concerning water-graphite contact angle and in our 
opinion they better reproduce the conditions represented in 
molecular models of such systems. Taking this into consid- 
eration, we decided to calculate new Lennard- Jones param- 
eters for water-carbon interactions, using the recent flexible 
SPC/Fw water model B45I1 . We basically followed the pro- 
cedure reported by Koumoutsakos and co-workers, although 
with some slight modifications which we will describe in the 
next section. Moreover, along with spherical droplets, we also 
simulated cylindrical filaments [46] (see Figure [1} of equiva- 
lent size and compared contact angle measurements for these 
two different systems. Defining a correlation between results 
for spherical droplets and cylindrical filaments has two impor- 
tant outcomes: first of all, it opens the possibility of avoiding 
time-consuming simulations of spherical droplets in further 
studies of static wetting properties of water on graphitic sur- 
faces, as already pointed out by other authors II26I1 . Secondly, 
by analyzing the different influence of droplet and filament 
size and shape on contact angle, it helps in understanding 
the molecular mechanisms determining nanoscale behavior of 
water on a graphitic surface. 



II. METHODOLOGY 



In order to identify the most suitable value for £co in the 
light of the recent experimental findings reported above, we 
followed the procedure adopted by Koumoutsakos and co- 
workers, i.e. we performed a series of MD simulations of 
water droplets of different radii on two graphene layers orga- 
nized as hexagonal graphite Il47ll and varied the value of £co 
in the range 0.1 —0.3 kJ/mol; then, we calculated the equilib- 
rium contact angle formed by the droplets on the surface and 
examined the relation between 9 and £co- The value of afco 
was kept fixed at afco = 3. 190 A iflill . We chose to use the in- 
creasingly widespread SPC/Fw water model ll45ll48ll49ll . for 
which the introduction of intramolecular degrees of freedom 
allows improved accuracy amongst the simple and common 
three point charge models. As previously said, both cylindri- 
cal filaments and spherical droplets were simulated, in order to 
understand how to relate contact angle measurements result- 
ing from these two different systems. Model size was varied 
in order to assess its influence on droplet or filament shape 
above the graphitic surface and to allow extrapolation of con- 
tact angles in the limit of very large base radius, i.e. for macro- 
scopic droplets. Finally, to confirm the behavior found using 
spherical droplets with sizes up to 8'000 molecules (3'213 for 
cylindrical filaments), we performed simulations of droplets 
composed of 16'000 water molecules (5' 117 for the equiva- 
lent filament) using £co = 0.20 kJ/mol. In the next sections, 
simulation and analysis details are reported. All computations 
were performed using Materials Studio 5.5 by Accelrys Soft- 
ware Inc. 
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this work 


Wu et al. (% difference) 


experimental (% difference) 


{roo) A 


1.0320 


1.03210 


0.1 


0.970 


6.4 


(0/hoh) 


107.47° 


107.69° 


0.2 


106° 


1.4 


p g/cm 3 


1.025 


1.012 


1.0 


0.997 


2.8 


O s lO^cmV 1 


2.24 


2.32 


3.4 


2.3 


2.6 



Table I: Bulk water properties: comparison and % difference with data from Wu et al. l46ll and experimental measurements (45 



III. SIMULATIONS 

Water model. We chose to base our study on the SPC/Fw 
water model. This model is strictly related to the widely 
used Simple Point Charge (SPC) water model introduced by 
Berendsen and co-workers 15011 . the main difference con- 
sisting of bond and angle flexibility within the single water 
molecule. Considering the SPC/Fw model (with roH P '"oh 2 
as the distance between the oxygen and each of the two hy- 
drogens and 0zhoh as the angle), the general interaction po- 
tentials can be expressed as 114511 : 



V" 



= y[('-OH 1 -^ H ) 2 

+ (>-OH 2 -roH) 2 ] + ^(eZHOH 



f/HOHJ 



v--TW[(3f)°-(a)1 + ^}. 



(3) 



(4) 



where V lntra and y lnter account for bonded and non-bonded in- 
teractions, Tq Hi and 0^ HOH are the equilibrium bond length 

and angle (1.012 A and 113.24° respectively), and k a 
are the spring constants (4431.534 kJ/mol/A 2 and 317.566 
kJ/mol/rad 2 respectively), ry is the distance between atoms 
i and j, Eij and ay, are the Lennard-Jones parameters for 
atom pair In the electrostatic term, qt is the partial 

charge on atom i (in units of elementary charge e). In the 
SPC/Fw model, no Lennard-Jones interactions are considered 
between oxygen and hydrogen, whereas oxygen-oxygen inter- 
action parameters are £oo = 0.650 kJ/mol and Obo = 3.165 
A. Partial charges on oxygen and hydrogen atoms are — 0.82e 
and +0.4 le respectively. 

In our simulations we deviated from the original work 
describing the SPC/Fw model and chose to use a differ- 
ent set of parameters for the evaluation of non-bond forces 
and for temperature and pressure control, more similar to 
those used in the works by Koumoutsakos and co-workers 
lfl8i l24l l25ll . Temperature and pressure control were achieved 
using Berendsen thermostat and barostat lf5lll while van der 
Waals and electrostatic forces were treated using group-based 
cutoffs. Cutoff length was set at r out = 10 A, with the poten- 
tial smoothed to zero using a switching function S(r) along a 
width of 1 A HI: 



S(r) 



2 

out 



Hr 2 



-2r 2 -3rl 



2 

out 



(5) 



where r; n = 9 A. 

Because of this choice of parameter set, it was necessary 
to check if our scheme could give reliable results in the cal- 
culation of bulk water properties. Therefore, we built a 4'000 
SPC/Fw water molecules box with periodic boundary condi- 
tions and we simulated its dynamical behavior in the NPT en- 
semble for 6.5 ns at a temperature of 298.15 K and a pressure 
of 1 atm, similar to what had been done in the original work 
J45|. We restricted ourselves to the calculation of basic prop- 
erties, i.e. average bond length (/"oh) and angle (0/hoh)> den- 
sity p, oxygen-oxygen radial distribution function g(roo) and 
self-diffusion constant D s . Results obtained for these analyses 
are reported in Table [I] along with comparisons with experi- 
mental measurements and results by Wu et al. l45ll . A graph 
of g(roo) (experimental and simulated) is reported in Figure 
[2] The slight differences between measured quantities and ref- 
erenced ones allowed us to regard our parameter set as being 
effective for the simulation of bulk liquid water at that specific 
temperature. Moreover, we could not spot any of the artifacts 
that could arise when using cutoffs for electrostatic terms in 
the simulation of water 1541 [55ll . Finally, as a final check for 
our settings, we calculated surface tension 7 for the SPC/Fw 
water model. Starting from an average density frame extracted 
from the last ns of the NPT trajectory, we performed an NVT 
simulation of a water slab by simply adding 50 A of vacuum 
space in the z direction. Dynamical behavior of the system 
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Figure 2: (Color online) Comparison between experimental and sim- 
ulation data for oxygen-oxygen radial distribution function. Experi- 
mental data are taken from Sorenson et al. 15311. 
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was again simulated for 6.5 ns at a temperature of 298.15 K, 
using a Berendsen thermostat and the previous cutoff scheme. 
Recording of pressure components P xx , P yy and P zz every 1 ps 
during the last 5 ns of the NVT simulation allowed us to cal- 
culate surface tension 7 by using the simple formula lf56l - l58tl : 

y= y (^-^), (6) 

where L z is the box length in the z direction. Surface tension 7 
estimated value was 70.8 ±1.9 mN/m, differing by only 1 .7% 
from the experimental value of 72 mN/m. 

Simulation of water on graphitic surface. Once our parame- 
ter set was proven to hold reliable results for water bulk prop- 
erties and surface tension, we built models of water spher- 
ical droplets and cylindrical filaments, starting from an aver- 
age density frame selected from the NPT simulation described 
above. From this state, we extracted several hemispherical and 
hemicylindrical shapes, whose characteristics are reported in 
Table [TT] For the largest models, we simply replicated the 
4'000 molecule box in the three directions of space, in order 
to have enough molecules for our structures. 

Hemispherical and hemicylindrical shapes, containing a 
number of water molecules ranging from 250 to 16'000, were 
then positioned just above (approximately 2 — 3 A) the model 
of two graphene sheets arranged as hexagonal graphite, using 
periodic boundary conditions. We used only two graphene 
sheets as long as additional ones would have been beyond the 
cutoff distance for any water molecule on the surface; more- 
over, graphene sheet atoms were kept fixed during all sim- 
ulations, consistently to what had been reported in previous 
works 1 18]. Side lengths of the graphene sheets were varied 
according to the droplet shape and size. For spherical droplets, 
square planes were used, with sizes approximately ranging 
from 100 x 100 to 160 x 160 A 2 . For cylindrical droplets, 
we used rectangular planes, with a fixed thickness of 25.6 A 
and a length ranging from 100 to 160 A. The length of the 
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Figure 3: (Color online) Average density map for a spherical droplet 
comprising 8'000 water molecules and £co = 0.2 kJ/mol. Color code 
based on density expressed as number of molecules per A 3 . Similar 
representations are obtained for the other droplets, both spherical and 
cylindrical. 



simulation box in the z direction was kept equal to the longest 
edge in the xy plane. Box boundaries were distant enough to 
prevent any interaction between periodic images for any of the 
simulated models. 

All of the systems described above were simulated using 
different values of £co- Simulation procedure was similar to 
the one used by Werder et al. 118I1 and consisted of three dif- 
ferent steps. Firstly, a simple energy minimization was per- 
formed; secondly, for equilibration purposes, each system was 
simulated in the NVT ensemble for 0.5 ns at a temperature of 
298.15 K, using a timestep of 1 fs. A Berendsen thermostat 
was used for temperature control and non-bond interactions 
were treated with the scheme described above, i.e. with a cut- 
off distance r out = 10 A and a smoothed potential calculated 
using Equation|5] Finally, a 0.5 ns simulation was performed 
in the NVE ensemble, again with a timestep of 1 fs and the 
same non-bond interaction scheme. 

Analysis. In order to calculate the contact angle of spher- 
ical droplets and cylindrical filaments of water on graphene 
sheets, we first had to define an equilibrium profile. For this 
purpose, we used two variations of the analysis procedure re- 
ported by different authors Il8ll36ll37ll . For spherical droplets, 
the method basically consists of dividing the volume occu- 
pied by water molecules by using a cylindrical binning, reg- 
ular along the z axis and in the angle 9. By imposing that 
the length of the n-th radial bin is r„ = ro\/« + 1, all elements 
have the same volume dV, ro being the length of the first ra- 
dial bin; in this work, we chose dV = 0.009 A 3 . Reliability 
of this choice was proven by considering different yet reason- 
able values for dV as well as for length and angle segmenta- 
tion, which led to almost identical final results for the contact 
angle measurement. It is worth noting that some care should 
be taken in order to avoid finite size effects. For every frame 
of the NVE trajectory (frames were taken every 0.5 ps), the 
number of atoms falling within a given volume element were 
counted and averaged over both the total number of frames 
and the angle 0, thus obtaining an equilibrium density map 
for the droplet (Figure |3). Water bulk density pb was calcu- 
lated in the inner part of the droplet density map, i.e. for those 
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Figure 4: (Color online) Equilibrium profile calculated with cylin- 
drical binning procedure (black points) and circular fit (red line) for 
a spherical droplet comprising 8'000 water molecules and £qo = 0.2 
kJ/mol. Tangent line determining contact angle is shown in green. 
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spherical droplets 


cylindrical filaments 


radius r\ A 


No. water molecules 


radius r\ A 


No. water molecules 


14.9 


250 


14.9 


322 


18.6 


500 


18.6 


490 


23.8 


l'OOO 


23.8 


737 


30.1 


2'000 


30.1 


1'263 


38.5 


4'000 


38.5 


2'068 


48.1 


8'000 


48.1 


3'213 


60.7 


16' 000 


60.7 


5' 117 



Table II: Initial radius r\ and number of water molecules of droplet and filament models. 
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l/r„ [1/A] l/U [1/A] 

Figure 5: Left: Cosine of the contact angle 6 as a function of 1 /r ca for spherical droplets. Right: Cosine of the contact angle as a function 
of 1// C a for cylindrical filaments. Error bars indicate standard deviation calculated by averaging over measurements taken on 5 blocks of 100 
ps each during the NVE production run. 



volume elements located from h a /4 to 3h a /4 along the z axis 
and from to rf/4 along the radial axis (h a and rf being the ap- 
proximate droplet height and radius obtained from the density 
map). For every layer in the z axis, then, we detected those 
volume elements for which density first exceeded pb /2, start- 
ing our check from the vapor phase. In this way, we obtained 
equilibrium profiles which were subsequently approximated 
by a circular fit (Figure |4). Conversely to what had been done 
by other authors II 18ll36l 13711 . we decided to include all points 
in the fitting in order not to introduce any arbitrary cutoff for 
the definition of the circular profile. We also decided to use 
a circular fit in order to keep the whole procedure as simple 
as possible, even though other investigators preferred to em- 
ploy more sophisticated fits B5711 . Here we restrict ourselves 
to a circular fit because experimental measurements are gen- 
erally carried out under the assumption of a spherical shape of 
droplets. 

For cylindrical filaments, we adopted an analogous proce- 
dure except for the fact that in this case averaging was done 
over the whole thickness of the periodic box instead than over 
the angle 9. In any case, we kept dV = 0.009 A 3 also for 
these systems, in order to allow a comparison as rigorous as 
possible between data obtained from cylindrical filament and 
spherical droplet simulations. 



IV. RESULTS AND DISCUSSION 

In Figure [5] cosines of contact angles calculated for 
spherical droplets are plotted against contact line curvature 
1 /r ca . A similar plot for cylindrical filaments is also reported, 
where Z ca is intended as the inverse of half the length of the 
chord determined on a vertical section by the points at the in- 
tersection between the circular fit and the plane (see also Fig- 
ure [TJ. By plotting cosines of contact angles against these 
two quantities, we intend to provide a term of comparison be- 
tween droplets and filaments with the same £co and initial 



£co kj/mol 


0oc 


A0o„ 


spheres 


filaments 


0.10 


157.0° ±6.9° 


147.9° ±5.2° 


9.1°±12.1° 


0.15 


132.6° ±4.5° 


140.8° ±6.4° 


8.2° ± 10.9° 


0.20 


123.2° ±5.0° 


128.9° ±6.8° 


5.7° ±11.8° 


0.25 


118.5°±3.3° 


116.6° ±3.4° 


1.9° ±6.7° 


0.30 


108.0° ±4.0° 


108.7° ±5.7° 


0.7° ±9.7° 



Table III: Extrapolated values of 6L for different values of eco- 
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♦ E=0.30 [kJ/mol] 
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Figure 6: Contact angle as a function of rf for spherical droplets, Left, and cylindrical filaments, Right. Error bars indicate standard deviation 
calculated by averaging over measurements taken on 5 blocks of 100 ps each during the NVE production run. 
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Figure 7: 0oo as a function of £co for spherical droplets (cf. Table 

UnJ. 



radius. 

From Equation I2 we would expect a monotonous trend in 
the curves of Figure[5] but for spherical droplets two different 
regimes seem to be present. For larger droplets, contact angle 
slightly decrease with increasing r ca and it is possible to de- 
termine an average positive line tension of (5.060 ± 1.061) • 
10~ 12 N, in line with results reported in literature lfl8l [35I1 . 
Conversely, for very small droplets (fi < 30 A), contact an- 
gle seems to decrease with decreasing r ca . In any case, if we 
consider only those droplets with r; > 30 A, we can easily 
extrapolate the macroscopic contact angle 9„ for the different 
£rjo values we considered (see TableHTTb. In order to match the 
experimental value of 127°, an appropriate estimation for this 
parameter for SPC/Fw water should be around 0.2 kJ/mol (see 
Table IllTl and Figure [7]). This is significantly lower than what 
is commonly used in oxygen-carbon Lennard- Jones potentials 
for modeling the interaction between water and graphitic sur- 
faces lfl7H29ll . assuming water is modeled with a set of pa- 
rameters similar to the one described above (i.e. three point 
charge models with a Lennard- Jones potential between oxy- 




Figure 8: (Color online) Schematic representation of the section 
of two spherical droplets/cylindrical filaments with different radii 
r\ and (section is intended perpendicular to the graphitic sub- 
strate and to the cylinder axis). The dotted line represents the in- 
teraction range of a water molecule close to contact line. The two 
droplets/filaments can not possibly present the same value for the 
contact angle (as is depicted), not only because of line tension ef- 
fects related to contact line curvature, but also because the resultant 
force acting on the molecule in the smaller droplet/filament lacks the 
contribute of the water molecules in the shaded area. 



gen atoms). 

For cylindrical filaments, conversely, the value of contact 
angle should be independent from their size, as long as contact 
line curvature is zero for this geometry and therefore no con- 
tribute should be present from the term related to line tension 
in Equation [2] Nonetheless, we can see a slight increase in 
with Z C a (see Figure|5]), similarly to what is noticed for smaller 
droplets. Table [Till reports contact angles for filaments if this 
trend is extrapolated to the macroscopic limit; as it would be 
expected, results show good agreement with the correspond- 
ing contact angles calculated for spherical droplets. Interst- 
ingly, if we plot contact angle versus rf (i.e. the final radius 
of the sphere/cylinder) for both filaments and droplets (see 
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Figure |6j, it is possible to see that 9 in small droplets it is 
more sensitive to size variations via rf than its counterpart in 
cylindrical filaments. The reason for this unexpected behavior 
in our MD simulations for both filaments and small droplets 
(i.e. decreasing 9 for decreasing droplet/filament size, hence 
contact angle values not consistent with modified Young's 
law) could reside in the fact that, in these systems, molecules 
near the contact line experience reduced cohesive forces be- 
cause of the higher curvature of the liquid/vapour dividing 
surface and this results in an increase in the adhesive com- 
ponent of the total force acting on them. This effect can be 
easily understood from the schematic representation of Fig- 
ure [8] where the interaction range of the molecule is depicted 
as a circle of radius r ou t- This reduction in cohesive forces 
has an opposite effect to the one related to a positive line ten- 
sion (which for spherical droplets should in principle increase 
contact angle with decreasing droplet size) and it seems to be 
dominant over the one related to contact line curvature 1 /r ca 
for spherical droplets with r\ < 30 A. A possible validation for 
this supposition comes from the fact that for cylindrical fila- 
ments, where no effect of line tension is expected, the increase 
in contact angle with increasing size is still present and it is 
less marked than for equivalent small spherical droplets, as it 
would be logical considering the difference of these surfaces 
near the contact line. Moreover, it should be considered that as 
the droplet/filament grows smaller, the number of molecules 
which belong to the most superficial layer of the droplet in- 
creases relatively to those in the bulk. This could again result 
in the reduction of cohesive forces for molecules near the con- 
tact line and in the decrease of the contact angle. One could 
argue that these results are merely an artifact caused by the 
choice of using group or atom based cutoffs, but additional 
simulations using larger cutoffs (12.5 and 15 A) for £co = 0.2 
kJ/mol gave similar results (data not shown), as already re- 
ported by other authors for atom based cutoffs H25I1 . even if 
in that case convergence of contact angle measurements was 
reached for slightly larger values of r ou t- Similarly, it has been 
proven that more accurate methods for the treatment of long 
range interactions in this kind of MD simulations do not have 
a significant influence on the value of contact angle ll24ll . Fi- 
nally, it is worth noting that both spherical droplets and cylin- 
drical filaments, simulations with the lowest value of £co cor- 
respond to more scattered results and deviations from the be- 
havior noted for all other values of the interaction parameter 
(see Figures [5j ID |7]), indicating that when adhesive forces are 
extremely low it becomes more difficult to unambiguously de- 
fine and measure contact angle for the simulated systems. 



V. CONCLUSIONS 

We have run MD simulations of spherical droplets and 
cylindrical filaments of SPC/Fw water on two sheets of 
graphene arranged as hexagonal graphite and calculated the 
relative equilibrium contact angles along the contact line. Us- 
ing system parameters similar to those previously adopted in 
the literature for the simulation of spherical droplets, results 
indicate that, in order to recover the recently measured macro- 
scopic contact angle of 127° between water and graphene, 
the value of £co in the Lennard-Jones potential used to de- 
scribe non bonded interactions between oxygen and carbon 
should be around 0.2 kJ/mol, much lower than the one that 
is commonly used in MD studies. Comparison between con- 
tact angle measurements for spherical droplets and cylindri- 
cal filaments seems to indicate that these two systems lead 
to comparable results in the macroscopic limit, but both sys- 
tems show that contact angle varies with drop size in a way 
which is not consistent with modified Young's equation. As 
already suggested by other authors [59], it seems that, beyond 
line tension correction, other contributions should be con- 
sidered in modifying Young's equation for extremely small 
droplets/filaments; we believe these contributions could be re- 
lated to the dividing surface curvature and to the relative in- 
crease in the number of molecules belonging to the most su- 
perficial layer with respect to the bulk phase. In any case, 
it has to be noticed that these results have been obtained in a 
very specific framework and, to a certain extent, are dependant 
on the choice of the parameters used in the simulation. More- 
over, due to the challenges presented by measurements of con- 
tact angle for such small systems, it would be very difficult to 
seek an experimental confirmation to the reported behavior. 
Therefore, we believe that further work in MD simulation of 
this sort is necessary, in order to sistematically investigate the 
influence of all different parameters and thus to help in the 
understanding of wetting phenomena at the nanoscale. 
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